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Abstract 

We evaluate Shannon entropy for the position and momentum eigenstates of some conditionally exactly 
solvable potentials which are isospectral to harmonic oscillator and whose solutions are given in terms 
of exceptional orthogonal polynomials. The Bialynicki-Birula-Mycielski (BBM) inequality has also been 
tested for a number of states. 
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1 Introduction 



In recent years there have been a growing interest in studying information theoretic measures for quantum 
mechanical systems. In particular entropic uncertainty relations (EUR) which serve as alternatives to 
Heisenberg uncertainty relation have been examined by a number of authors pQ. Among the various 
measures of information theoretic entropy, the Shannon entropy plays a particularly important role and 
Shannon entropy of the probability distribution is used as a measure of uncertainty. An EUR relating 
position and momentum was obtained by Beckner, Bialynicki-Birula and Mycieslki (BBM) [5] and it is 
given by 

Spos + S mom > D(l + log tt) (1) 
where D denotes the spatial dimension and the position and momentum entropies are defined by 

S P os = -! \^{x)\ 2 log\4>(x)\ 2 d D x, S mom = - J \^{p)\ 2 log\^(p)\ 2 d D p (2) 

where ip(x) is a normalized coordinate space wave function and ip(p) is the corresponding Fourier trans- 
form. Apart from their intrinsic interest the EUR ([1} or the information entropies ^ have been used 
in different contexts e.g, to study squeezing [3], localization and fractional revivals [4j [5] etc. However 
it is not always easy to exactly evaluate Shannon entropies ([2]) or the EUR (p}. They have been ob- 
tained exactly only for a few low lying states of the harmonic oscillator [TJ [5] , Poschl- Teller [7] [5] , Morse 
potential |7j (9], Coulomb potential etc [10] and for other states especially those with large quantum num- 
bers numerical results have been obtained [10] . In a related development position space entropy for some 
potential isospectral to the Poschl- Teller potential has also been computed 11 . Also information en- 
tropies of several classical orthogonal polynomials e.g, Laguerre, Hermite, Gegenbauer etc., which provide 
solutions of most of the standard solvable potentials have been calculated [12] . However there are many 
exactly solvable potentials, especially the isospectral partners of the standard ones whose solutions can 
not be expressed in terms of classical orthogonal polynomials and as far as information theoretic measures 
are concerned not much is known about these potentials. 

Very recently a new class of orthogonal polynomials, called the exceptional orthogonal polynomials 
[13] . have been studied by a number of authors [14]. These polynomials are different from the classical 
orthogonal polynomials and have rather distinct properties. Interestingly exceptional orthogonal poly- 
nomials multiplied by square root of the associated weights appear as solutions of a number of exactly 
solvable potentials. Here we shall examine a class of conditionally exactly solvable potentials [15] which 
are isospectral to the harmonic oscillator potential |16[ I17j and whose solutions are given in terms of 
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exceptional orthogonal polynomials related to the generalized Laguerre and Hermite polynomials |18j . 
These potentials consist of a harmonic oscillator term and other terms of a non polynomial type and 
may be viewed as a deformation of the harmonic oscillator 19 . In view of the fact that the isospectral 
partner potentials have rather different type of cigenfunctions it is of interest to evaluate their position 
and momentum space entropies and compare them with the same of their partners. In this paper we 
shall consider the models of ref [18] and examine them from the point of view of EUR (H} . 

2 Exceptional orthogonal polynomials associated with Condi- 
tionally exactly solvable potentials 

It may be noted that a supersymmetric quantum mechanical system consists of a pair of Hamiltonians 
(H ± ) of the form (in units of % = m = 1) [20] 

H ± =A±A* = -~+V±(x), A ± = ±=(±± + W{x)}, V ± {x) = \(W\x)±W\x)) (3) 

Thus zero energy states of the Hamiltonians H± can be obtained from the solutions of the equations 
A T ipy = and they are given by = e ± / WP(a),te . If either of ^ is normalizable supersymmetry is 
unbroken while if neither is normalizable then supersymmetry is said to be broken. However both of ip^ 
can not be normalizable at the same time (at least in the examples we shall consider) . So if supersymmetry 
is unbroken then the above Hamiltonians are isospectral except the zero energy state (which is assumed 
to belong to In this case the relationship between the energies and the eigenfunctions of the 

Hamiltonians H ± are given by 

E = 0, E- +1 =E+>0 (4) 
i>o=Ne~f w ^ dr , ^+ = _^A+^ +1 , C+t = 4= A~rP+ (5) 



e- +1 y/K 

Thus all the states are doubly degenerate except the zero energy ground state. On the other hand 
when supersymmetry is broken we have 

E+ = E->0, i;+ = -±= A+tl>-, il)- = -^=A-^+ (6) 
V E n v E n 

Thus in this case all the states including the ground state are doubly degenerate. We shall now examine 

the case of broken supersymmetry and consider the following superpotential : 

/ + 1 n'(r) 

W(r)=r + — + ^4, 0<r<oo, ;>0 (7) 
r u(r) 
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where u(r) is given by 



u(r) = xFt 



1 -e 



(8) 



2 ' 1 2' 

It is seen that u(r) will be a polynomial for any positive odd integral value of e. It can be checked that 
neither of tp^ is normalizable and consequently supersymmetry is broken. Now we choose the simplest 



non trivial value i.e, e = 3 so that 



Then it can be shown that [I7J EH] 



u(r) 



V+(r) 



E+=E-= 2n + 21 + 5, = 



2r 2 + 21 + 3 
21 + 2, 



1(1 + 1) , 7 

_: i J- 7 J 

2r 2 + + 2 



2(n!) 



r(n + Z + f) 



r z+1 iL +i (r 2 )e- r2 / 2 



, n = 0, 



(9) 

(10) 
(11) 



y (r) = ^ + ( f + 1 )( f + 2 ) + ^ 

" W 2 2r 2 (2r 2 + 2? + 3) 



2r 



2(1 + 1) 



4r 



2r 2 + 21 + 3 



g 2 y-^ - ! - ^ 



4iL +i (r 2 ) + 2(2r 2 + 21 + 3)L l n l (r 2 ) 



(12) 



(13) 



y (2n + 21 + 5)r(n + i + |) (2r 2 + 2Z + 3) 
There are several points to be observed. First to be noted is that the solutions (fT3| can be expressed in 
terms of the following exceptional orthogonal polynomials and the weight function in the following way 



4>n ( r ) = N ny /w{r) p n (r) 
Pn(r) 



w(r) 



(21 + 3) 
2Z + 3 



4L^ +5 (r 2 ) + 2(2r 2 + 21 + 3)l'+ 2 (r 2 ) 



(14) 



2r 2 + 21 + Z, 

Secondly V+ (r) is the effective potential for the radial oscillator while V- (r) is the conditionally exactly 
solvable isospectral partner which is non shape invariant. 



3 Information entropy for the potential VL(r) 

It may be noted that so far we have considered a system on the half line and I is just a parameter. 
However if we wish to view our system as a three dimensional system with a radially symmetric potential 
then I has to be treated as angular momentum. It may be observed that the radial part of V+ (r) is that 

1 If a potential is exactly solvable only if some of the couplings in the potential assume specific values, then such a 
potential is termed as conditionally exactly solvable | 15| . 
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of a standard radial oscillator while that of V- (r) actually depends on I and consequently the degeneracy 
pattern for a fixed value of I = L is given by L+l = E~ +1 L . Nevertheless it will be seen that the 
operators A + or A~ may be used to simplify calculation of the position entropy of the (— ) sector. We 
would like to note that V+(r) represents the effective potential for the radial harmonic oscillator. On the 
other hand the potential V— (r) is isospectral partner of V+ (r) . However V_ (r) is of anharmonic type and 
it is not shape invariant. Here our objective is to evaluate Shannon entropy for the isospectral partner 
V-(r) and examine their differences with those of V+(r). 

The potential V- (r) in ([12")) and the wave function (| 13[) are the effective potential and the reduced wave 
function. The actual spherically symmetric potential is the same as (|12[) except the angular momentum 
term while the complete wave functions are given by 



= V (2n + 2/ + 5)I> + / + l) WT2TT3) ^ ^ + ^ + % + ^ ^ \ 

(15) 

where YJ™ 1 (0, </>) denote the angular part of the wave functions. The position entropy can be readily 
evaluated using (|15l) . To evaluate the momentum entropy it is necessary to obtain the Fourier transform 
of (|15[) . In this context it may be noted that momentum space wave functions for the harmonic oscillator 
or the Coulomb potential can be obtained in a relatively simple manner |22j . However for other potentials, 
particularly the isospectral partners of the standard potentials, the momentum space wave functions can 
not be obtained in the same way. Here our approach would be to express the momentum space wave 
functions of V- (r) partly in terms of the momentum space wave functions of the radial oscillator which 
are already known [TOl [22] . To this end we note that the Fourier transform of the wave function (|15[) is 
given by 22 

fc+l,tn(P) = ^+l(P) i H , l(e.*) ( 16 ) 



where 



91+1(P) = Jl I™ 3l+i(prWn,i dr, ji(z) = J^J l+i (z), (17) 



Ji+l{z) being the Bessel function (of the first kind) of order (I + |) [21] ■ Now using the intertwining 
relation ^ inside the integral (|17[) we obtain 



, -C+i) / 32 n l 

V2n + 2i + s < I+1 , m (p) + — —\l nr{n + l+ s ) ^?i(e,*) Hp) 



(18) 
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where t/£ j+1 m (p) is the Fourier transform of V^~z+i m( r ) QUI 122] : 



2n! 



/ T(n + i + §) 

The function 7(p) is the solution of the differential equation 



p l+l e -^/2 £ ^3/ V) y l!pi(e)$) 



dp 2 



g±iW±2 + ( , + ') 



f 



(19) 



= (20) 



and is given by 



J(p) = VP 



J i+3 / 2 HV^+372p) (ci + ( l J^' 2 ^%- f2 / 2 ll +3/2 (-zV^+372^)^ +5/2 (< + jft 2 ) + Li +§ (i 2 ))dt) 



2+3/21 



A-iy/T+mP) (ga - ( ^ 3/!> ^ e- t2/2 J ;+ 3 /2 (-tVT+372t)t' +5 / 2 (4 + _|(t 2 ) + (t 2 ))dt 

(21) 

where F;(z) denotes Bessel function of the second kind [21 . The constants C\ t i depend on the quantum 
numbers n, I and have to be determined from the boundary conditions 



(22) 



It may be pointed out that the expression (|18j) for the momentum eigenstates is an exact one. Now 
using (fT5|) and (fT8|) we shall compute the Shannon entropies ([2]) and examine the EUR ([lj for different 
values of the quantum number n and the results are given in Tables 1 and 2. We would like to note that 
in Tables 1 and 2 we have considered those eigenfunctions which correspond to the same energy. From 
Table 1 it is seen that the BBM inequality is always satisfied as it should be and the sum of entropies 
increases with n. From Table 2 we find that except for n — 0, for the same energy the sum of the entropies 
for the isospectral partner (fT"2)) is less than the same for the original potential. 



3.1 Information entropy for isospectral partner of linear harmonic oscillator 

Here we take the superpotential to be jTSJ [T7] 



W(x) = x + 



Ax 



1 + 2x 2 ' 

Then the partner potentials and their solutions are given by 



oo < x < oo 



1 



2™n!0r 



-x 2 /2 



H n (x), E+ = n + 3, n = 0,l--- 



(23) 



(24) 



(» 



x 2 4 16a; 2 3 

{X) l + 2x 2 + (l + 2x 2 ) 2 + 2 : 




e -x 2 /2 

~ (1 + 2x2)' E o =0 



(25) 



-x 2 /2 



[(1 + 2x 2 )H n+1 (x) + AxH n {x)} , E- +1 =n + 3 



v / 2"+ 1 n!(n + 3)0F (1 + 2a; 2 ) 
In this case the polynomials and the weight function are given by 

Po(x) = 1, p n (x) = (1 + 2x 2 )H n+1 (x) +4xH n (x), oj(x) 



(26) 



(l + 2ir 2 ) 2 

It may be noted that V+(x) in (|25p is a linear harmonic oscillator (with a shifted energy scale) while its 
partner is a conditionally exactly solvable potential. A major difference with the previous example is that 
in the present case V— (x) has a zero energy ground state and consequently supersymmetry is unbroken. 
Thus the ground state of V-(x) is a singlet while the excited states are degenerate. Also because of 
the relations (|5|) it is clear that the ground state V'o" is annihilated by both A + and A~ and has to be 
considered separately. Following the procedure of the last section it can be shown that 



A (p) 



(ire) 



1/4 



2V2 



2cosh(^=) - erf[\ - JL] 



er f[l 



V2 J 



where 



y/2(n + 3)2 n n!0F 

V3+(p) = 



1 



v^F2«+i(n + l)! $+ +1 (p) + U(p) 
e-v 2 l*H n {p) 



(27) 
(28) 

V*^~ ~ ny " (29) 

is the momentum space wave functions of the potential (|24l) while I(p) is the solution of the differential 
equation 

^T(n\ 1 ,-n-l p -p 2 /2 

(30) 



d 2 I(v) 1 i n - 1 e~ p 1 2 

—Jl - -I(p) + ==— [2nH n ^(p) - P H n {p)\ = 
dp 2 2 ^2 n+2 n\y/TT 



and is given by 

— p 

I(p) — e 73 

,-n+l 



C 2 + e 
■p 



s/2p 



,-n+l 



2 V /2"+ 1 n! % /7T Jl 



(2nH n -i(q) - qH n (q))dq 



2 y /2 n + 1 n\ifK Ji 



±q(V2-q) 



(2nH n -i(q) - qH n (q))dq 



(31) 
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As in the previous example the constants C'1.2 depend on the quantum number n and have to be deter- 
mined for each n from the condition 

lim ^" +1 (p) = (32) 

p— >±oo 

Using (f24]) -(f2"8 |) we have calculated position and momentum entropies for a few levels and the results are 
given in Table 3. From Table 3 it is seen that in contrast to V-(x), the ground state of the isospectral 
partner does not saturate the BBM inequality. Also for states corresponding to the same energy position 
and momentum entropies of the isospectral partner is more than the same of the linear oscillator. 

4 Conclusion 

In this paper we have computed the Shannon entropies for the position and momentum eigenstates for 
a first few levels of two conditionally exactly solvable potentials whose solutions are given in terms of 
exceptional orthogonal polynomials. It has been found that for linear oscillator the entropies of the 
eigenstates of the isospectral potential are more than those of the linear oscillator. However in the radial 
problem the situation is similar for the initial levels although the angular momentum of the eigenfunctions 
differ by one unit. 

Here we have considered a certain type of potential and we feel it would be of interest to examine other 
potentials whose solutions are given in terms of exceptional orthogonal polynomials of other types e.g, 
exceptional Jacobi polynomials |14j . It may also be noted that we have obtained the results for the some 
initial states and it is desirable to be able to evalaute entropies for any value of the quantum number(s). 
However computation of entropies for large quantum numbers is a formidable job. In view of this we 
feel it would be interesting to study asymptotic behaviour (with respect to the quantum number) of the 
entropies of exceptional orthogonal polynomials. Finally we would like to mention that it would be worth 
investigating some other properties like the Fisher information [23], spreading length [24 etc. of this new 
class of polynomials. 
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n 


Spos ( n ) 




Spos( n ) + S mom (n) 


3(1 + log(n)) 





3.361 


3.646 


6.917 


6.434 


1 


4.015 


4.199 


8.214 


6.434 


2 


4.568 


4.628 


9.196 


6.434 


3 


4.822 


4.954 


9.776 


6.434 



Table 1. Information Entropies and their sum for the 
eigenstates of V- (r) in (|12|) for I = m = 0. 



n 


Spos ( n ) 


Spos( n ) 









3.217 


3.361 


3.217 


3.646 


1 


4.151 


4.015 


4.151 


4.199 


2 


4.709 


4.568 


4.709 


4.628 


3 


5.109 


4.822 


5.109 


4.954 



Table 2. Information Entropies for the 
eigenstates of the potentials in (|10p and (112[) 

for I = m = 0. 
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Spos( n 1) 




Spos ( n 1 ) 






Spos{ n ) 


1 + log{-K) 











0.479 


1.679 


2.158 


2.144 


1 


1.072 


1.072 


2.144 


1.261 


1.607 


2.868 


2.144 


2 


1.343 


1.343 


2.686 


1.425 


1.578 


3.003 


2.144 


3 


1.499 


1.499 


2.998 


1.578 


1.748 


3.326 


2.144 



Table 3. Information Entropies and their sum for the eigenstates of the potentials 
in (|24j) . 



9 



References 

[1] R.J. Yafiez, W. Van Assche and J.S. Dehesa, Phys.Rev. A50, 3065 (1994) 
W. Van Assche, R. J. Yafiez and J. S. Dehesa, J.Math.Phys. 36, 4106 (1995) 

A. I. Aptekarev, J. S. Dehesa and R. J. Yafiez, J. Math. Phys. 35, 4423 (1994). 

[2] W. Beckner, Ann.Math. 102, 159 (1975) 

I. Bialynicki-Birula and J. Mycielski, Comm.Math.Phys. 44, 129 (1975) 
I. Bialynicki-Birula and L. Rudnicki, preprint arXive:1001.4668. 

[3] A. Orlowski, Phys.Rev. A56, 2545 (1997). 

[4] R. Atre, A. Kumar, C.N. Kumar and P. Panigrahi, Phys.Rev. A69, 052107 (2004). 
[5] E. Romera and F. dc los Santos, Phys.Rev.Lett. 99, 263601 (2007). 
[6] V. Majernik and T. Opatrny, J.Phys. A29, (1996) 2187. 

[7] J.S Dehesa, A. Martmez-Finkelshtein and V.N. Sorokin, Mol.Phys. 104, 613 (2006). 
[8] M.W. Coffey, Can.J.Phys. 85, (2007) 733. 

[9] E. Aydiner, C. Orta and R. Sever, Int.J.Mod.Phys. B22, 231 (2008). 
[10] J.S. Dehesa, W. Van Assche and R.J. Yafiez, Mcthod.Appl.Analysis. 4, 91 (1997). 
[11] A. Kumar, Ind.J.Pure.Appl.Phys. 43, 958 (2005). 

[12] J. S. Dehesa, R. J. Yafiez, A. I. Aptekarev and V. Buyarov, J.Math.Phys. 39, 3050 (1998) 
V.S. Buyarov, J.S. Dehesa and A. Martinez-Finkelshtein, J. Approx. Theory. 99, (1999) 153 
V.S. Buyarov, P Lopez-Artes, A Martmez-Finkelshtein and W. Van Assche, J.Phys. A33, 6549 
(2000). 

[13] D. Gomez-Ullate, N. Kamran and R. Milson, J. Math.Analysis.App. 359, 352 (2009). 

[14] C. Qucsnc, J. Phys. A41, (2008) 392001; SIGMA. 5, 084 (2009) 

S. Odakc and R. Sasaki, Phys.Lett. B679, (2009) 414; ibid, B682, 130 (2009); ibid, B684, 173 
(2010); J.Math.Phys. 51, 053513 (2010) 
T. Tanaka, J.Math.Phys. 51, 032101 (2010) 

B. Midya and B. Roy, Phys.Lett. A373, 4117 (2009) 



10 



B. Bagchi, C. Quesne and R. Roychoudhury, Pramana J.Phys. 73, 337 (2009) 

C. L. Ho, S. Odake and R. Sasaki, preprint, larXiv:0912.5447l 

D. Gomez-Ullate, N. Kamran and R. Milson, preprint arXiv:1002.2666 

[15] A. de Souza Dutra, Phys.Rev. A47, R2435 (1993) 

R. Dutt, A. Khare and Y.P Varshni, J.Phys. A28, L107 (1995) 

R. Roychoudhury, P. Roy, M. Znojil and G. Lcvai, J.Math.Phys. 42, 1996 (2001). 

[16] G. Junker and P. Roy, Phys.Lett. A232, 155 (1997). 

[17] G. Junker and P. Roy, Ann.Phys. 270, 155 (1998); Phys.Atom.Nucl. 61, 1736 (1998). 

[18] D. Dutta and P. Roy, J.Math.Phys. 51, 042101 (2010). 

[19] D. Gomez-Ullate, N. Kamran and R. Milson, J.Phys. A37, 1789 (2004). 

[20] F. Cooper, A. Khare and U. Sukhatme, Supersymmetry in Quantum Mechanics (World Scientific, 
2001) 

G. Junker, Supersymmetric Methods in Quantum and Statistical Physics, Springer- Verlag, 1996. 
[21] I.S. Gradshteyn and IM. Ryzhik, Table of Integrals, Series and Products, Academic Press, 1994. 
[22] S. Fliigge, Practical Quantum Mechanics, Springer- Verlag, 1994. 

[23] R.J. Yahez, P. Sanchez-Moreno, A. Zarzo and J.S. Dehesa, J.Math.Phys. 49, 082104 (2008) 
J.S. Dehesa, S. Lopez-Rosa, B. Olmos and R.J. Yahez, J.Math.Phys. 47, 052104 (2006). 

[24] A. Guerrero, P. Sanchez-Moreno and J.S. Dehesa, J.Phys. A43, 305203 (2010) 

P. Sanchez-Moreno, J.S. Dehesa, D. Manzano and R.J. Yahez, J.Compt.App.Math. 233, 2136 (2010). 



11 



